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Abstract 

We study the constraints on dark energy equation of state u x and the deceleration parameter q from 
the recent observational data including Hubble data and the cosmic microwave background (CMB) radiation 
by using a model-independent deceleration parameter q(z) = 1/2 — a/(l + z) b and dark energy equation of 
state u x = ojo + u)iz/(l + z) in the scope of anisotropic bianchi type I space-time. For the cases of Hubble 
dataset, CMB data, and their combination, our results indicate that the constraints on transition redshift z, 
are 0.62lJ;gg, 0.34tg;J|, and O.bCpQ^g respectively. 
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1 Introduction 

According to the recent cosmic observations our universe is experiencing an accelerating phase of expansion 
at the present time m-0- Since the ordinary barionic matter (energy) can only produces an attractive force, 
there should necessarily an exotic form of energy with negative pressure called “dark energy” (DE) be exist 
in order to drive the observed cosmic acceleration m-m • Moreover, based on the recent observations we live 
in a nearly spatially flat Universe composed of approximately 4% baryonic matter, 22% dark matter and 74% 
dark energy. A natural candidate for dark energy is cosmological constant A, which has the equation of state 
wa = pa/ PA = —1. However, although cosmological constant can explain the present cosmic acceleration, it 
would encounter some serious theoretical problems, such as the fine-tuning and the coincidence problems. An¬ 
other possible candidates for dark energy are the dynamic dark energy models provided by scalar fields such 
as quintessence — | > w > — i i-m, phantom ( u x < —1) |T5] , quintom (uA < —|) [15], k-essence [T7, IB], 
Chaplygin gas as well as generalized Chaplygin gas models nans], and etc. 

We also can study dark energy through an almost model-independent way. In this approach we parameterize 
the dark energy EoS parameter by giving the concrete form of the equation of state of dark energy directly, 
such as w(z) = ujo + uqz [S], uj(z) = wo + wi [52, :53], etc. Studies of dark energy as a function of redshift, 
by using this method, show that though the cosmological constant is not ruled out in lcr region, the current 
constraints favors a dynamical dark energy. However, as noted in ref [53] the rapid changed of EoS parameter 
i.e <C 1 is ruled out. Moreover, it is a well established fact that the universe is accelerated expansion at 
present (dark energy dominated) and decelerated expansion in the past (dark matter dominated). Therefore, the 
deceleration parameter (DP) should not be a constant but time (or redshift) variable which is positive in the past 
and changes into negative at z* ~ 0.5 [24]-[26[. Considering the idea of variable DP, the parameterized deceler¬ 
ated parameter is present in almost model independent way by giving a concrete form of decelerated parameters. 

The high symmetry involved in FLRW models requires a very high degree of fine tuning of initial conditions 
which is extraordinary improbable and hence FLRW models are infinitely improbable in the space of all possible 
cosmologies m- Although there are an increasing interest in the study of DE in the scop of anisotropic 
space-times, there is not a reference paper including the present observational constraints on the observable 
parameters such as q , w, z*, etc. Therefore, still authors compare their results by those obtained on the bases 
of the standard cosmological models. Motivated the situation discussed above, in this letter instead of FLRW 
metric, an anisotropic space-time namely Bianchi type I metric is used in order to obtain more general results. 
Up to our knowledge, such kind of calculations has not been done yet. First we derive the general form of 
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the EoS parameter in Bianchi type I space-time and then we use the Hubble dataset, CMB data, and their 
combination to put the observational constraints on the model parameters. I is worth to mention that because 
of technical difficulties in minimizing \ 2 the results obtained in this work are not so general and more work is 
needed to extend these results for example by minimizing Xsneia- 


2 Dark Energy Equation of State 

We consider the Bianchi type I space-time in the orthogonal form as 

ds 2 = — dt 2 + A 2 (t)dx 2 + B 2 (t)dy 2 + C 2 (t)dz 2 , 
where A(t),B(t) and C(t) are functions of time only. 

The Einstein’s field equations ( in gravitational units 8nG = c = 1) read as 


(1) 


m - - Rg) = T 


( m)i 


T, 


{X)i 


( 2 ) 


where 2y m ^' and are the energy momentum tensors of barotropic matter and dark energy, respectively. 

These are given by 


Tj n ' i = diag[—p m ,p m ,p m ,p m ], 
= diag[—1, w m ,w m , w m ]p” 


(3) 


and 


Tj = diag[—p ,p ,p ,p A ], 


* ~ X :P X dp X ] 

X , X , .X-I X 


= diag[-l,w ,w ,w A ]p , (4) 

where p^ m ^ and are, respectively the energy density and pressure of the perfect fluid component or ordinary 
baryonic matter while w m = p m /p m is its EoS parameter. Similarly, p x and p x are, respectively the energy 
density and pressure of the DE component while w A = p x /p x is the corresponding EoS parameter. We assume 
the four velocity vector u l = (1,0, 0,0) satisfying u l Uj = — 1. 

In a co-moving coordinate system (u l = <5q), Einstein’s field equations (0 with 0 and 0 for B-I metric 0 
subsequently lead to the following system of equations: 


^ + ^ + ^ = - Ujmpm -“ XpX ' 


A C AC m m y x 
A + C + AC = ~“ p - UJp ' 
A B Ab m m xx 
A + B + AB=^ P -“ P ' 


AB AC BC 
AB + AC + BC~ r+P ' 


_m | 


(5) 

( 6 ) 

(7) 

( 8 ) 


If we consider a = (ABC) 3 as the average scale factor of Bianchi type I model, then the generalized mean 
Hubble’s parameter H defines as 


TT a 1 A B C 
H ~a~3\A + B + C 


(9) 


The Bianchi identity Gf- = 0 leads to Tf- = 0. Therefore, the continuity equation for dark energy and 


baryonic matter can be written as 

p m + 3H(1 + uj m )p m +p x + 3H(1 + co x )p x = 0. 
Solving eqs. 00 one can find [28] 

A(t) = a±a exp(b\ J ci~ 3 dt), 


( 10 ) 

( 11 ) 
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( 12 ) 


and 

where 

Using eqs. m-m in eqs. (El¬ 
and 


B(t) = d 2 d exp(b 2 j d 3 dt ), 
C(t) = a 3 d exp(b 3 J d~ 3 dt), 


(13) 


aia 2 a 3 = 1, b 3 + b 2 + b 3 = 0. 

i we can write the analogue of the Friedmann equation as 


P 


Ka 


-6 


= -g(P + 3p). 


(14) 


(15) 


Here p = p m + p de , p = p m + p de and K = bib 2 + b\b 3 + b 2 b 3 . Note that K denotes the deviation from isotropy 
e.g. K = 0 represents flat FLRW universe. Thus, when the universe is sufficiently large, almost at the present 
time, the space-time © behaves like a flat FLRW universe. 

Moreover, if there is no interaction between dark energy and cold dark matter(CDM) with w m = 0, then one 
can write the conservation equation m for CDM and dark energy separately as separately as 


and 

Ea. dTfil) leads to 

Using eqs. m, mi in eqs. © and 

and 


p x + 3H(l+co x )p x = 0, 
p m + 3Hp m = 0. 


p m = p™d~ 3 . 


we obtain the energy density and pressure of dark fluid as 
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P X = 3 


H 

Ho 


- 3 Kd~ 6 - p™a 


m 3 


p x = -2- - 
d 


H 

Ho 


— Ld 


(16) 

(17) 

(18) 

(19) 

( 20 ) 


respectively. Here, L = b% + b 3 + b 2 b 3 is a positive constant (Note that K + L = 0). Therefore, the equation of 
state parameter (EoS) of DE in it’s general form is given by 


u x = P 


, X _ 2q — 1 — La~ e (y^) 


-2 


y™n -3 


3 + 3La -6^ '-3Hg 
where q = is the deceleration parameter and Hq 1 is the current value of matter density. 


( 21 ) 


3 Experimental Tests 

In this section, by using x 2 method we investigate the constraints on the parameters ui x , q , and z*. To do so we 
utilize the recent observational data including Hubble parameter and the cosmic microwave background (CMB) 
radiation. Note that as mentioned by Kumatsu et al [29] the baryonic oscillation (BAO) distance ratio is not 
applicable in non-FLRW based models. 
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First of all, from eq. m we obtain the expression for the Hubble rate, of the form 


2 = & L(1 + 3cu x )(l + z)~ 6 

Hi (2g-l)±3uW(H™(l±z) 3 -l)' 

In our study, following Xu L X et al. SUL the deceleration parameter is considered to be 


q(z) = 


1 

2 


a 

(1 + z) b ' 


(22) 


(23) 


where a and b are constants determined by the recent observational constraints. From this equation it is obvious 
that for z 1, q — > i which is corresponding to matter dominated era whereas for z = 0, the current value of 
deceleration parameter is obtained as qo = ^ — a. In view of eq. (1231) . the Hubble parameter is written in the 
form 

H(z) = H 0 ( 1 ± z)^exp [a((l + z)~ b - 1 )/b] , (24) 

We also consider the following simple expression for the equation of state 

U) X =UJq +idi ~—■—, (25) 

I ± z 

where to o and uq are constants. Using eqs. (E51) and (051) one can easily re-write eq. E21) as a function of z only. 


As far as we are dealing with the experimental H{z) test, it has been suggested by Jimenez et. al. (31, to use 
the quantity 4| which is called differential age. In this method, one can directly analyze the H(z) data without 
passing through the luminosity distance as we usually do in the case of Supernovae la. Using the differential 
ages of passively-evolving galaxies from the Gemini Deep Deep Survey (GDDS)[32j and archival data [531 [54] . 
first, a set of nine values of H{z) in the range of 0 < z < 1.8 has been reported [55]. Table 1 shows the most 
recent data set includes 14 values. 


We minimize the following reduced % 2 [36] to constrain the parameters oj x and q for fixed Hq at 71 and 

(26) 


flfi 1 = 0.27. 




2—1 


'Hub 


0 *) 


where H obs are the values of Table 1 and v = 13. 


Table 1: The cosmological data at lcr error for H(z) expressed in s 1 MPc 1 Km. 


z 

0.00 

0.10 

0.17 

0.24 

0.27 

0.40 

0.43 

0.48 

0.88 

0.90 

1.30 

1.43 

1.53 

1.75 

H(z) 

72 

69 

83 

79.69 

77 

95 

86.45 

97 

90 

117 

168 

177 

140 

220 

la error 

±8 

±12 

±8 

±4.61 

±14 

±17 

±5.96 

±60 

±40 

±23 

±17 

±18 

±14 

±40 


To constraint our models we have also used the CMB. In this case, by using the shift parameter R one can 
constraints the model parameters by minimizing 


9 \R th - R obs V 

XCMB ~ 


(0.019) 5 


Here R obs = 1.725 ±0.018 


and R th is given by 

R th = (n™)i 


dz 

m 


(27) 


(28) 


where zcmb = 1091.3. 

The constraints on model parameters can be obtained by minimizing XHub^~XcMB- First we rise eq. (PU1) in eqs. 
(EHD and m to find constraints on a, b 1 and q , then we apply eqs. (E51) and m on eq. (1^1) using the results 
from previous calculations, by ftting the observational Hubble data, CMB, and their combination respectively, 
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Table 2: The best fit parameters with ler error. 


observational data 

"V Z . 
a. mm 

a 

b 

w 0 

W 1 

Qo 

z* 

Hubble data 

9.02 

u - y '-0.56 

i q'7+U.t>/' 

± -°' -0.25 

-1-08-1.O2 

0.83!^ 

-o.47i8:S 

0 62 -1 " 1 ' 45 

u - uz -0.56 

CMB data 

161.21 

1.43t°;g 

2 in + 0 - 89 
1.60 

-1.42±8;®i 

l-25±S;i| 

-0.93li s 46 

0 34 + 0.13 
u -°^-0.06 

Combination 

168.34 

1-26 

1 Q 6 +0 ' 50 

l.yu_i 89 

-o.98t<|;l| 

0-43±g;ll 

-o.76i8:H 

0 60 +0 ' 20 
u - uu -o.io 





(a) Hubble (b) CMB (c) H+CMB 

Figure 1: Plots of the parameters a and b with ler confidence level according to: (a) Hubble, (b) CMB, and (c) 
H + CMB combination data 


we obtained the values of minimum x , best fit parameters a , b wo and wi; transition redshifts z* , and current 
deceleration parameter qq. Our results are listed in Table 2. 

The ler contours of parameters a and 6 , the variation of the decelaration parameter q , and the ler contours of 
parameters uq and w 0 are depicted in figures 1, 2, and 3 respectively. All these Figures show that the combined 
constraint is stricter than independent constraints given by CMB dataset or the Hubble data. 

The best fit transition redshift z* for the case of combined constraint at the ler error is 0.601 qiq. It is worth to 
mention that the constraints to the deceleration parameter q and the dark energy equation of state are looser in 
ler error by using Hubble data than CMB data, which may be attributed to a fewer data points from Hubble data. 

The variations of the DP parameter q and the EoS parameter ui x are plotted in fig. 4. Here we have used 
the best fit values of a, b , wq, and uj± from Table 2. 





(a) Hubble (b) CMB (c) H+CMB 

Figure 2: Evolutions of the decelerated parameter q versus red shift z, which are constrained according to: (a) 
Hubble, (b) CMB, and (c) H + CMB combination data 
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(a) Hubble 


(b) CMB 


(c) H+CMB 


Figure 3: Constraints on the dark energy equation of state oj x = wo + wiy^ obtained by: (a) Hubble, (b) 
CMB, and (c) H + CMB combination data. A cross sing in figures shows cosmological constant corresponds to 
w x = -1. 



z 


Figure 4: The plot of DP parameter q and H x versus redshift (z) for a = 1.26, b = 1.96, to o = —0.98, and 
wq = 0.43. These values are taken from the combined data in Table 2. 
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